Magnetic field generation in relativistic shocks 
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Abstract. We present an analytical estimate for the magnetic field strength generated by the Weibel 
instability in ultra-relativistic shocks in a hydrogen plasma. We find that the Weibel instability is, 
by itself, not capable of converting the kinetic energy of protons penetrating the shock front into 
magnetic field energy. Other (nonlinear) processes must determine the magnetic field strength in the 
wake of the shock. 

The following article has been submitted to 2003 GRB Conference Proceedings. After it is pub- 
lished, it will be found at http : / /proceedings . aip . org/. 



INTRODUCTION 

The fireball model for Gamma-ray Bursts (GRBs) explains GRB afterglows as syn- 
chrotron radiation coming from the external shocks that are formed when a relativis- 
tically expanding fireball (or jet) interacts with surrounding gas [1]. The spectra and 
luminosity of the observed afterglows indicate a magnetic field strength in the radiat- 
ing material of about 10% of the equipartition field strength [2]: B 2 /Sn ~ £ with e the 
total post-shock energy density. Such a magnetic field is much stronger than what is ex- 
pected from simple shock physics: for instance, passive compression of a dynamically 
unimportant magnetic field in a relativistic shock yields a post-shock field satisfying 
(for example, see [3]) B 2 /Sne ~ (va/c) 2 , with va the Alfven speed ahead of the shock, 
which usually satisfies va <C c. 

It may be possible that much stronger magnetic fields are generated in the shock 
transition itself through a Weibel-like instability. This low-frequency electromagnetic 
beam-instability can develop at the point where the shocked and unshocked plasma 
penetrate each other. It converts the kinetic energy of the penetrating particle beams 
into thermal motions and unordered magnetic fields [4]. For external shocks propagating 
into a hydrogen plasma incoming protons carry most of the kinetic energy in the shock 
frame. The instability must convert a significant fraction of this energy into magnetic 
energy in order to reach the required field strength [5]. Here we investigate if this is 
indeed possible. 



SIMPLE COLLISIONLESS SHOCK MODEL 

Within the shock transition, where the fireball material encounters the surrounding 
plasma, the mixing of the shocked (relativistically hot) and the unshocked (cold) plasma 
produces a plasma with a very anisotropic velocity distribution. As in non-relativistic 



collisionless shocks, it is quite likely that a significant fraction of the incoming ions 
is reflected at the shock transition by a large-scale electrostatic or magnetic field. The 
resulting situation is known to be unstable [4]: in the so-called Weibel instability the 
penetrating and reflected particles will bunch together, causing electric currents that 
induce a magnetic field in the plasma, which in turn causes the particles to bunch even 
more. If the generated magnetic fields become sufficiently strong they will trap the beam 
particles, eventually saturating the instability. 

Recent numerical simulations [6], [7], [8] show that the electric currents generated 
by the Weibel instability merge with each other in the wake of the unstable region. This 
separate process has a longer time-scale, and forms structures on a larger length-scale 
than the Weibel instability, which occurs mostly on a scale of the order of the plasma 
skin depth cj (0p e , with (0p e the electron plasma frequency (see below). 

In this paper we will concentrate on the proton-driven Weibel instability. 



THE WEIBEL INSTABILITY AND ITS SATURATION 

Because of the small electron mass, the electron-driven Weibel instability evolves very 
rapidly [5], with the magnetic field strength growing as e at with a « (Opb where 

C0pb = ^J^Tlehi^Jm^ is the beam particle plasma frequency based on the proper beam 
density n 0h . When the electron-instability has saturated, the electron velocity distribu- 
tion will be close to isotropic. The electrons then form a relativistically hot background 
plasma in which the much slower proton-driven Weibel instability develops because the 
proton velocity distribution is still very anisotropic (see, for example, figure 6 of [7]). 
We will investigate the behavior of the protons in this situation. 

The main features of the resulting instability can be reproduced by investigating a 
water-bag proton velocity distribution. We assume two counter- streaming proton beams 
moving along the jc-direction, with a small velocity spread in the z-direction to model 
thermal motions. The proton momentum p is then distributed as: 

F(p) = -r^- [§ (p x - Pxo) + S(Px + Pxo)]S(p y ) [0 (p z - p z q) - ®(p z + pa)], (1) 

where n v is the total proton density, p x Q is the beam momentum, p z o is the maximum 
momentum in the perpendicular direction, 8(x) is the Dirac delta function and 0(x) 
is the unit step function. We consider the evolution of a wave perturbation with wave 
vector k = ke z and frequency ft). The growth rate of the proton instability is calculated 
in the usual manner by looking for wave solutions of the linearized equations of motion 
and Maxwell's equations (for example [9]). The resulting dispersion relation links the 
complex wave frequency ft) to the wave number k. The Weibel instability obeys a 
dispersion equation of the form (using the notation of [9]): 

(O 2 -k 2 c 2 +C xxz = 0, (2) 

where C xxz contains contributions from both the electrons and the protons. Since the 
electrons are relativistically hot, their contribution is C xxz ^ = —C0 2 e , with (Dp e the effec- 
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FIGURE 1. Left: growth rate as a function of wave number for a shock with Lorentz factor 1000 and 
thermal velocity spread v,o = 0.001 . The solid line is the result for an electron-proton plasma. The dashed 
line is the result for pure pair plasmas. Right: magnetic field strength as a function of wave number for a 
shock with the same parameters and n p = 2 cm -3 . 



tive electron plasma frequency which equals 

C0pe = \j l 4ne 2 n e /m e h , (3) 

where h = (e + P) e /n e m e c 2 is the electron enthalpy divided by the rest mass energy 
density, which parameterizes the relativistic mass correction in a relativistically hot 
plasma. 

The proton contribution is [9]: 

c„ = «a ^ (T) + J_ *J_ > , (4) 



-po y Yb j Yb0(O 2_ k 2 v 2 



zO 



where C0p[ = ^Ane 2 n v /m v is the (non-relativistic) proton plasma frequency based on 
the lab-frame beam density, m p is the proton rest mass, u\ = pi/ (m p c), fbo = (1 +"x0 + 

M Jo) 1/2 ' v ' = u i c /7b0 and 



i ln fi±V£, (5) 



With these results the dispersion equation (2) becomes a biquadratic equation for 0) 
that has one positive imaginary solution giving the growth rate o = Im(ft)) of the 
unstable mode (figure 1, left). 

One expects the instability to saturate when the quiver motion of the beam particles in 
the wave reaches an amplitude Az such that kAz ~ 1 . An expression for the corresponding 
magnetic field strength is derived in [4] and reads in our notation: 

5 S „ = ^^. (6) 

Using the dispersion relation between o and k we can then find 5 sat as a function of k 
(figure 1, right). 



A straightforward calculation shows that for small v z q the peak of 5 sa t lies at wave 
number fc peak given by 

^peakC 2 - + (7) 

with a corresponding growth rate 

The proton-driven Weibel instability occurs for a relatively small range in wavelength 
compared to the electron-positron case (figure 1): modes with wavelength longer than 
the electron skin depth (k < c/a>pe) are inhibited by the response of the background 
electrons to the proton perturbations. The location where the growth rate levels off 
corresponds with the location where the saturation magnetic field peaks (k = fc pea k)- If the 
background electrons had not been as responsive, this location might have been at much 
lower wave number, and since the wave number is in the denominator of expression (6) 
for 5 S at, this would have resulted in a much higher magnetic field. However, we find 
that the electrons set the location of the peak (fc pea k — Wpe/ c) and that the proton Weibel 
instability can only produce slightly stronger magnetic fields than the electron instability 
(figure 1, right), despite the larger kinetic energy of the protons. 



CONCLUSION 

The Weibel instability as it was modeled here is not efficient at converting the kinetic en- 
ergy of the protons into magnetic fields. However, numerical simulations of collisionless 
shocks in electron-proton plasmas ([6], [7], [8]) do show efficient production of mag- 
netic fields. This can probably be attributed to merging of the electric currents produced 
by the Weibel instability, after the instability that was considered here has stopped. To 
verify this, the current merging processes will need to be investigated further. 

The properties of the magnetic fields generated by these processes will be important 
for determining the radiation that may be produced in these shocks and that we see in 
the form of gamma-ray burst afterglows. 
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